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ABSTRACT 


The  purpose  of  this  report  is  to  explicitly  derive  two  sets  of  six-degree- 
of  freedom  equations  of  motion  for  a  symmetric  missile. 

One  set  is  based  upon  a  coordinate  system  that  is  rigidly  attached  to 
the  missile  (body-fixed  system),  while  in  the  second  set  (fixed-plane  sys¬ 
tem)  a  coordinate  system  with  one  axis  constrained  to  lie  in  a  given  plane 
is  employed  to  derive  the  equations  of  motion. 

Both  sets  of  equations  assume  the  earth  to  be  spherical,  include  the 
effect  of  the  earth's  rotation,  and  consider  variable  wind.  In  addition,  for 
the  body-fixed  system,  discussion  of  initial  conditions  and  singularities 
of  the  differential  equations  is  presented. 
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TABLE  OF  SYMBOLS 


Coordinate  systems 


Arbitrary  vector  quantities 


Angular  velocities 


V  Missile  velocity,  independent  of  wind  velocity 

-# 

V*  Wind  velocity 

Vf  Missile  velocity  relative  to  the  wind 


if 

-♦ 

XL 

-* 

J 


Ix»  !y*  lz 

0,  ^ 

R 


Ar 

-• 

Y 


Summation  of  forces  acting  on  missile 
Summation  of  torques  acting  on  missile 
Total  angular  momentum  of  missile 

Moments  of  inertia  about  the  three  axes  of  the  missile 
Euler  angles 

Position  of  missile  center  of  gravity  relative  to  reference  co¬ 
ordinate  system 

Increment  of  time 

Unit  vector  indicating  direction  of  gravitational  force 
Vector  indicating  altitude  of  missile 
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Mean  radius  of  the  earth 


Gravitational  factor 
Mass  of  rocket 
Density  of  air 
Diameter  of  missile 

Distance  from  aose  of  missile  to  center  of  pressure,  expressed  in  calibers 

Distance  from  nose  of  missile  to  msgnus  center  of  pressure,  expressed 
in  calibers 

Distance  from  noae  of  miaaile  to  missile  center  of  gravity,  expressed 
in  calibers 

Thrust  malalignment  distances 

Thrust  of  missile 
Axial  drag  coefficient 
Normal  force  coefficient 
Magnus  force  coefficient 
Yaw  damping  moment  coefficient 

Roiling  moment  coefficients 

Proportionality  factor  of  thrust  applied  to  jet  torque 
Angle  between  yaw  plane  and  >M  **** 

Angles  giving  position  of  thrust  relative  to  missile  coordinate 
system 
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A0  Longitude  of  missile  launch  position 

B°  Latitude  of  missile  launch  position 

H°  Missile  angle  of  declination  relative  to  geographic  system 

G°  Angle  measure  d  from  positive  k"c  axis  relating  lateral  missile  posi¬ 

tion  in  geographic  system 

_  Denotes  time  derivative 

dt 

Nntotlen 

A  general  vector  V  will  be  represented  in  a  particular  coordinate  system  by 
(V),  where  i  denotes  the  referencing  coordinate  system. 

A  component  of  a  vector  will  be  denoted  by  two  subscripts,  the  first 
giving  the  component,  the  second  denoting  the  coordinate  system  used. 

-* 

Here  X,  Y,  and  Z  represent  components  along  the  i,  j,  k  axes  respectively.  For 

-♦ 

eiampie  (Vw  }y£  is  the  component  of  the  wind  velocity  along  the  j  £  axis. 

Finally  ^  denotea  differentiation  relative  to  the  ith  coordinate  system, 
dt 
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INTRODUCTION 


This  report  is  part  of  a  confirming  program  to  give  Picatinny  Arsenal  a 
complete  capability  in  the  flight  simulation  of  all  types  of  projectiles  and 
missiles,  whether  ballistic  or  rocket-boosted,  guided  or  unguided.  This 
capability  is  important  to  the  development  of  both  conventional  and  spe¬ 
cial  weapons.  It  provides  necessary  information  for  aerodynamic  design, 
range  calculations,  error  analyses,  fuzing  systems,  and  complete  weapon 
systems  evaluations.  The  mathematical  model  defined  in  this  report  forms 
the  basis  for  a  corresponding  digital  computer  program  which  has  been 
developed  for  the  IBM  709/1401  system. 

Contained  herein  are  two  derivations  of  the  six-degree*of-freedora  equa¬ 
tions  of  motion  of  a  missile.  The  two  derivations  differ  by  the  fact  that, 
in  one  case,  forces  and  moments  acting  on  a  missile  are  referred  to  a  co¬ 
ordinate  system  that  is  rigidly  attached  to  a  missile  (body-fired  system), 
while,  for  the  second  set,  the  forces  and  moments  are  referred  to  a  co¬ 
ordinate  system  .that  has  one  axis  constrained  to  lie  in  a  given  plane 
(fixed-plane  system).  It  is  believed  that  the  first  set  is  more  appropriate 
for  asymmetric  missiles  with  the  missile  coordinate  system  chosen  to 
coincide  with  the  principal  axes  of  the  body.  This  system,  however,  seems 
to  have  a  disadvantage  for  high  spin  rates  (spin-stabilized  rockets)  in  that 
it  may  be  necessary  to  take  very  small  time  increments  to  obtain  an  accu¬ 
rate  trajectory  by  numerical  methods  (this  will  be  indicated  in  RESULTS 
AND  DISCUSSION).  For  this  case  the  second  set  of  equations  seems  more 
appropriate. 

It  should  be  noted  that  it  may  be  advantageous  to  choose  other  coordi¬ 
nate  systems  than  considered  here,  either  to  suit  a  particular  missile  or 
because  of  the  type  of  resultant  data  desired.  Stability  considerations  and 
terminal  effects  are  examples  for  which  specialized  results  are  required. 

Both  sets  of  equations  in  this  report  include  the  effect  of  the  earth’s 
rotation  about  its  axis,  gravitational  expressions  for  both  the  flat  earth 
(the  earth  is  assumed  to  be  a  plane,  valid  for  short  trajectories)  and  the 
spherical  earth  (the  earth  is  assumed  to  be  a  perfect  sphere),  variable  air 
density,  and  the  effect  of  wind.  It  is  also  well  to  state  specifically  that 
the  following  effects  are  not  included  in  the  derivations: 
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1.  Guidance  factors 


2.  The  motion  of  the  earth  along  its  orbit 

3.  Launcher  effects 

4.  Asymmetric  missiles 

5.  Stability  criteria 

These  represent  areas  of  extension  to  the  present  equations.  The  reader 
is  referred  to  References  1,  2  and  3,  which  treat  in  some  detail  some  of 
these  additional  factors. 
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To  introduce  the  reader  to  the  equations  of  motion,  a  brief  derivation  is 
presented  in  Equations  1  to  11.  The  definitions  of  the  various  symbols 
are  given  on  pages  4,  5,  and  6.  Equations  1  and  2  provide  the  foundation 
for  the  equations  of  motion.  Arrows  over  symbols  indicate  vector  quantities. 


Here  the  derivatives  are  taken  with  respect  to  a  fixed  coordinate  system 
(ij,  ij,  kj).  Now 
-# 

—  -  v  «  VXIij  +  Vyjii  +  VZjkj  -  VXyiM  +  VYMiM  +  vZUkU  (3) 
dt 

iM,  in.  k y  indicate  a  moving  coordinate  system. 

Equation  1  can  be  written  in  terms  of  the  moving  coordinates  as 


„4i. 

dt  L  dt  J 


where  _  indicates  a  derivative  relative  to  the  moving  coordinate  sys¬ 


tem.  In  addition. 


a>  *  <^XI‘I  +  ^YIJI  +  ^Zl^l  *  UYM‘M  +  WYM*M  +  u  Zkfc  M  (5) 


which  upon  expansion  yields 

ZF=  m  d  VxM  iM  +  +  ,.dVz-k  +  (*>YMVZM  ~  WZMVYM>‘M 

dt  M  dt  M  dt  ^ 


+  (wzmVXm  ~  <*>XMVZM>iM  +  <  WXMVYM  “  wYMVXM>kM 


in  terms  of  the  moving  coordinate  system. 
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Similarly  for  Equation  2 


2l  +(;*;> 

dt  dt 


(7) 


where  r.ow 


J  "  +  +  IZMwZMkM 


(8) 


which  upon  expansion  yields 


£.wu±3Su.ivU-l^HL 


IXM  lyM - Lfi  j*  *  *Z»4  - “kM+  wyMWzM^ZM  -IYm)'U 

dt  dt  dt 


♦  ^XM^Zh^XM  ~  ^ ZM^  M  +  ^XJ^YM^YM  _IXM>kM 


(9) 


also  in  terms  of  the  moving  coordinate  system. 

The  forces  and  moments  are: 

£F  -  (Fx  +  Fjj  gravitational  +  FX  chruat^M 

-♦ 

+  (^Y  normal  +  ^Y  magnua  +  ^Y  gravitational  +  ^Y  thmat^M 
+  <FZ  normal  +  FZ  magnua  +  FZ  gravitational  +  FZ  thruat^kM  (10) 

EL  =  (Lx  toll  damping  +  jet  torque  *  LX  apin  deceleration 

-♦ 

+  Lx  thruat  malalignment^M 

+  (Ly  restoring  +  Ly  yaw  damping  4  Ly  magnua 
+  Ly  thrust  malalignment^!  M 

+  (LZ  restoring  +  Lz  yaw  damping  +  Lz  magnua 

+  Lz  thrust  malalignment^  M  (11) 
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In  succeeding  sections  of  this  report,  these  equations  of  motion,  along 
with  all  the  necessary  supplementary  equations  for  their  solution,  will 
be  explicitly  derived.  In  particular,  Part  A  of  the  section  entitled  PRO¬ 
CEDURE  covers  principles  of  rotating  coordinate  systems  and  Euler 
angles.  In  Part  B,  these  principles  are  used  in  deriving  the  equations  of 
motion.  A  complete  tabulation  of  the  resulting  equations  is  presented  at 
the  end  of  the  PROCEDURE. 

PROCEDURE 

Several  coordinate  systems  are  used  in  deriving  the  equations  of  motioo. 
The  reasons  can  be  summarized  as  follows: 

1.  Forces  and  moments  acting  on  a  missile  are  commonly  resolved 
along  the  axes  of  a  "missile  coordinate  system,"  that  is,  a  coordinate 
system  with  a  conveniently  specified  orientation  relative  to  the  missile 
itself.  It  should  be  noted  that  this  coordinate  system  will  usually  travel 
in  some  manner  along  the  missile  trajectory. 

2.  Newton' s  laws  of  motion  are  strictly  valid  in  an  inertial  coordi¬ 
nate  system,  that  is,  a  coordinate  system  fixed  in  space.  Because  of  this, 
any  forces  and  moments  that  are  resolved  along  the  axes  of  a  moving  co¬ 
ordinate  system  must  eventually  be  referred  to  fixed  (or  inertial)  coordi¬ 
nates  to  correctly  relate  forces  and  accelerations. 

3.  It  is  often  convenient  to  refer  the  motion  of  a  missile  to  coordi¬ 
nates  that  are  neither  inertial  nor  dependent  upon  missile  orientation.  An 
example  is  the  case  when  it  is  desired  to  refer  the  motion  of  a  missile  to 
coordinates  fixed  on  the  earth’s  surface. 

For  these  reasons  the  basic  equations  governing  moving  coordinate  sys¬ 
tems  will  now  be  derived. 

A.  Mathematical  Preliminaries 

I.  Rotating  Coordinate  Systems 

Relations  that  exist  between  rotating  coordinate  systems  will  now 
be  derived.  Consider  first  two  non-rotating  right-handed  coordinate  systems 
(fA,  j a.  kA>  and  (fB.  jB,  kB>  as  in  Figure  1  <p  10).  A  vector  A  can  be 
represented  in  either  coordinate  system  as  follows: 
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PIV 


\V/s' 


*  <4 

Figure  1 


AA  -  AXA*A  +  aYAIA  +  AZAkA  (12) 

4  4  <4  «» 

AB  *  aXB‘B  +  AYB)B  +  AZB^B  (13) 

where,  since  a  vector  is  independent  of  the  coordinate  system  represent¬ 
ing  it, 

A  A  “  ab  *  A  (14> 

■4 

If  A  is  moving,  then,  relative  to  eac1-  coordinate  system 

,  d  *  A  *  d  A  v  a  d  Ay  !■*  d  A  7  . 

Va - - —  a+ - SlkA  (15) 


G. 

> 

>A 

> 

II 

d  Ay 

_  A  »  .  4- 

d  A  v  a 

IA  + 

dAZA 

dt 

A 

dt 

dt 

dt 

-♦ 

deAB  _ 

dAXB* 

d  Av  n-* 

—  “>6* 
dt 

dAZB 

dt 

j  B 

dt 

dt 

Since  the  time  rate  of  change  of  A  ( -  V)  is  another  vector  which  again  is 
independent  of  the  representing  coordinate  system, 


(17) 


Because  of  Equation  17  it  is  not  necessary  to  actually  specify 

or  to  which  coordinate  the  motion  of  A  is  calculated,  or 

dc  dt 
mathematically 


dAA 

dt 


dt 


dA 

dt 


(m 


This  is  not  t;  je  if  the  i„,  jR,  kB  coordinate  system  is  rotating  relative  to 

ddd  ^ 

iA,  j  A,  kA.  As  a  trivial  example,  it  is  possible  that  daA  •  0  while 

dt 

taL 

dt 


^  0,  for  the  case  of  a  vector  fixed  in  the  (iB,  jB,  kB)  system.  To 


calculate  the  rate  of  change  of  A  relative  to  the  fixed  coordinates,  but  in 

/ dA^B  \ 

terms  of  the  rotating  coordinates  I  -  I,  it  is  necessary  to  account 

'  dt  ' 

also  for  the  rotation  of  i  B,  jB,  kB,  namely, 


£*£§  =  ♦  .l£l£jB  «■  d^kB  +  AXBiAiB  +  Ay8  it!*  +  Azb (19) 

dt  dt  B  dt  B  XB  dt  YB  dt  ZB  dt 


To  determine  expressions  for  ,  ...  AfB  and  ^A^B  imagine  a 

dt  dt  dt 

vector  B  of  fixed  length  rotating  with  angular  velocity  w  as  shown  in 
Figure  2  (p  12). 
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Figure  2 

The  time  derivative  of  B  with  respect  to  any  fixed  system  is  as  follows: 

-+  +  +  - 
dB  .  Lim  B(t  ♦  At)  -  B(tl  „  AJL 
dt  Ar-»0  At  At-»0  At 

Referring  to  Figure  2,  the  magnitude  is 

i®|  .  *»B  «in  9  (20) 

At  I 

in  a  direction  perpendicular  to  the  plane  containing  &  and  t3.  This  vector, 
.  ...  .  ^  "• 
in  the  limit,  is  precisely  tu  ~  B. 

Consequently 

Lim  A®  *  il®  .  «  *  B  (21) 

At-*0  At  dt 


For  our  particular  case  B  can  represent  any  one  of  the  unit  vectors;  hence. 


d 


A*  B  _ 


dt 


=  cub  x  1  B> 


£^L®  =  CJBxjB;  ^  A.k  B.  =  fc>B  *  h 


dt 


B  *  *B 


(22) 


where  <uB  is  the  angular  velocity  of  ifi,  jB,  kB 
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Equation  19  becomes 


d  iAn  aA a  ** 

«  JLJL  +  axb(<d 


dt  dt 


AXB(  a>B  *  »B  )  ♦  A  yg<  wg  x  j  B  )  +  A zg(  <«>b  x  ^b)  (23) 


d«  A  &  "* 

where  _ Z— _L  denoc es  the  time  rate  of  change  of  A  relative  to  the  rotating 
dt  „ 

coordinate  system  (ifi,  jB,  kB).  This  derivative  shows  how  A  will  change 

in  time  to  an  observer  situated  on  the  rotating  coordinate  system. 
Rewriting  Equation  23  we  obtain 


dAAB  _  dsAfl 
dt  dt 

or  finally 


+  ( <uB  *  Aggig)  +  { «t*8  x  Aygj  g)  +  (<ug  x  A  2£kB) 


d  aAb  dBAB 


+  Wg  x  Xg 


This  is  the  basic  equation  governing  rotating  coordinate  systems. 

-+ 

To  determine  the  acceleration  of  A  relative  to  a  fixed  system,  again  in 
terms  of  the  rotating  coordinates,  we  have 


VAB 


(¥*-) 


Following  rules  just  established: 

-  Ig.  ( IaM  +(«bxV^b\ 

dt’  dt  '  dt  /  v  dt  I 


dB*AB  dB"B 


K  A  g  f  Cl)  g  x 


-»  -* 

+  (Og  x  -■  ■  ♦  *  (Wg  *  Ag) 

dt 
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Finally  we  have  the  result 


dA*AB 


■£  !tS.  +  d *  Ag  +  2tog  x  +  to g  x  (tog  *  Ag)  (26) 

dt*  dt  dt 


We  note  in  passing  that 


dA“n 

dt 


d  r.  to 


B  B 


dt 


+  ( tog  x  tog) 


or  since  to. 


toB  .  0 


dA6>B 

dt 


dBMB 

dt 


(27) 


(28) 


or  the  time  derivative  of  Z  can  be  taken  relative  to  either  coordinate  system. 
We  can  now  remove  the  restriction  that  the  two  origins  are  coincident. 

—a  -r 

Assume  that  the  (iB ,  }  B,  kg)  system  is  translating,  as  well  as  rotating, 
relative  to  the  (iA,  jA,  kA)  system.  We  have,  referring  to  Figure  3, 
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(29) 


«nd  finally 


dAC 

dt 


c 

-  A  ♦  B 

-* 

-*  <•* 

.  dAA  +  dAB 

dt 

dt  dt 

-* 

d*A 

dt 

dgB  -* 

♦  +  UU 

dt  B 

dAK. 


VA 


dt1  dt* 


B 

dt* 


'B  * 

2<w 


d„B 


-»  d. 


-*Ow  *»  .  *•  “•  Vi  _U>  _ 

B  *  -r~+  61b  X  <«B  *  B)  +  JLJ»  x  B 
dt  dt 


(30) 

(31) 


(32) 


This  discussion  is  concluded  by  a  theorem  that  will  be  of  considerable 
use  in  deriving  the  ballistic  equations  of  motion  when  more  than  two  co- 
ordinate  systems  are  involved. 

Given  a  primed  coordinate  system  rotating  with  angular  velocity  Z 

A 

with  respect  to  an  unprimed  coordinate  system,  and  a  starred  coordinate 
system  rotating  with  angular  velocity  <3  relative  to  the  primed  coordinate 

D 

system,  then  the  angular  velocity  of  the  starred  coordinate  system  relative 

to  the  unprimed  system  is  <wA  +  The  proof  is  as  follows:  The  velocity 
-+ 

of  any  vector  C,  fixed  in  the  starred  coordinate  system,  relative  to  the 
primed  coordinate  system  (denoted  by  d‘C/dt)  is 


d*  C 
dt 


Cl}  ^  X  C 


(33) 


and  the  velocity  of  C  relative  to  the  unprimed  coordinate  systems  becomes 


dC  d'C  - 


dt  dt 


.+  w1'<  C  *( x  C)  +  (  *  C)  =  ( Wg  +  )  *  C 


(34) 
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Before  deriving  the  equations  of  motion  a  second  mathematical  prelimi¬ 
nary  must  be  disposed  of,  namely,  the  Euler  angles  and  their  derivatives, 
it  is  presumed  that  the  reader  has  a  knowledge  of  the  elementary  operations 
of  matrices. 

II.  Euler  Angles  and  Their  Derivatives 

The  purpose  of  this  development  is  to  determine  what  relationships 
exist  between  the  axes  of  two  differently  oriented  coordinate  systems. 

These  relationships  are  important  since  we  may  know  the  components  of  a 
vector  (e.g.  velocity)  in  one  coordinate  system,  and  wish  to  know  the  com¬ 
ponents  of  the  same  vector  in  another  coordinate  system.  These  relations 
will  be  handled  by  what  are  known  as  Euler  angles,  although  one  may  use 
other  techniques  such  as  direction  cosines.  To  fix  ideas,  three  arbitrary 
rotations  of  a  (right-handed)  coordinate  system  will  be  performed  about 
selected  axes.  This  is  sufficient  to  orient  a  given  coordinate  system  into 
any  other  desired  (right-handed)  position. 

-4 

In  particular,  consider  a  fixed  coordinate  system  iA,  jA,  kA,  and  a  sec- 

ond  coordinate  system  (i,  j,  k)  initially  coincident  with  it,  and  whose  final 
position  will  be  denoted  by  ijj,  ja*  hg.  Intermediate  positions  of  i,  j,  k, 
will  be  denoted  by  a  sequence  of  primes,  the  number  of  which  denotes  the 
number  of  rotations  already  performed,  i  ,  .  .  .  ,  k  of  course,  all  being 
unit  vectors. 

Now  first  rotate  i,  j,  k  about  kA  by  an  angle  of  magnitude  as  shown  m 
Figure  4,  where  each  axis  of  the  triad  i,  j,  k  becomes  i*»  j*»  k*,  respectively. 


Ja 
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By  trigonometry,  the  projections  of  the  primed  axis  on  the  “A"  subscripted 
or  fixed  axis  can  be  determined.  The  result  of  this  computation  can  be  con* 
venientiy  written  in  matrix  form  as  follows: 


-♦ 

i' 

o 

o 

_ 1 

-b 

*A 

-♦ 

j* 

aa 

-S t/t  C0  0 

-♦ 

»A 

L  v  . 

.  0  0  1  . 

(35) 


It  is  important  to  note  that  this  matrix  equation  (and  all  similar  to  this)  is 
a  shorthand  notation  for  writing  three  equations  at  once.  Each  equation 
expresses  one  of  the  unit  vectors,  slong  the  coordinate  axis  of  one  co- 
ordinate  system  in  terms  of  the  other  (rotsted)  coordinates. 

For  the  sske  of  brevity  we  have  written 

SA  for  sin  A  and 
CA  for  cos  A 


where  A  is  any  angle. 

Similarly,  the  projections  of  the  unprimed  coordinates  on  the  primed  co¬ 
ordinates  can  be  written  as  the  inverse  (also  the  transpose)  of  the  shove 
coefficient  matrix,  namely, 


r  -  i 

“ 

lA 

C4>  -S0  o 

*' 

-b 

>A 

m 

S0  0 

i' 

«• 

-* 

.  kA  . 

0  ol 

k' 
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If  next  a  rotation  of  magnitude  6  about  j*  is  performed  we  obtain 


The  double  primes  are  used  to  denote  the  new  position  of  the  primed  co¬ 
ordinates. 
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Finally  we  perform  a  rotation  about  i"  of  angle  4,  obtaining 


and 


1  0  0 

i" 

*♦ 

ij» 

mt 

o  c<t>  S4> 

-* 

i" 

-4 

■*  j 

.  ^8 

o  -S4C4 

L  k"  J 

- 

* 

“ 

i#f 

1 

0  0 

*4 

j" 

0 

c<4  -s<p 

T 

>B 

*4 

kM 

0 

k 

s 4  c<t> 

-  *B  • 

(39) 


(40) 


Now  the  B-subscript  denotes  the  new  and  final  position  of  the  double 
primed  coordinates.  Three  sequential  rotations  of  magnitude  if>,  6,  and  4, 
respectively,  have  now  been  completed. 
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If  one  successively  multiplies  these  matrices,  the  following  is  obtained: 


-♦ 

*B 

cect 

cestf 

SB 

1b 

S«d>C0  -  C<£SySr 

sds<j>st  *  c$Ctp 

S4C6 

.  kB  - 

.  S6C<f>C<p  + 

sec -  'tyC'P 

c 4>ce . 

»A 

-♦ 

»A 


(41) 


and  the  inverse 


-  ^  ■ 

‘A 

cec# 

ses<j>ci!i  -  c<ps4 

«• 

iA 

m 

C.9 

SBS^Stp  +  C4>Cxp 

-* 

kA 

V 

.  SB 

S 4CB 

1 


c<t>c$ 


f? 


B 

*B 


L^B 


(42) 


To  obtain  the  significance  of  these  equations,  consider  a  vector  R 
written  in  terms  of  each  coordinate  system: 


R  -  (R)xA*A  +  (R)ya1a  +  (R)ZA^A  -  (R)XB»B  ■*  (R)YBiB  +  WzB^B  <42a> 


Assume  that  the  components  (R)xB’  (R)  Y B»  a0^  (^)zBare  knowo»  am*  *e 

wish  to  obtain  values  for  (R)xa*  (^)yA*  aoc*  ^ZA*  ”^°  °^ta^n  these  vaJ- 
-*  •«*  -¥ 

ues  we  write  iB,  jB,  and  k  g  in  terms  of  the  iA,  jA,  k  #  coordinates  and 
equate  like  components.  For  example,  from  Equation  41  we  have 


iB  *  CdC0iA  ♦  C5S^jA  -  Sfic  A 

■*  —  *  ■* 
j  B  -  (S0&J&C0  -  C«£Sy!0i  A  +  (S0S0S0  +  C4C0)jA  *■  S^CdkA  (42b) 

-*  - 

kB  -(S0C0C0  +S«iS0)iA  +  (S*:tfS0-S0C0>jA  +  CtfC0kA 

Substituting  these  values  in  Equation  42a  and  equating  the  coefficient^ 
i  on  the  right  hand  side  of  the  equation  to  (R)XA*  {h*  coefficient  of  iA 

on  the  left  side  of  the  equation,  and  the  coefficients  of  j  A  and  kA  to 
(R)ya  and  (R)za  respectively  yields: 

(R)xa  “  CBC<j(R)xb  +  (St&iiCit  -  C&Sv)  <r)yb  +  (SOCvCti  ♦  s<isd-)  (R)Za 

(R)  Y  A  =  COS^(R)  x  B  +  (S0Sd»Sti  +  C£C0)  (R)v  B  «■  (SOCpSii  -  StfCii)  (R)z  g 
(R)  ZA  =  -S^R)XB  *  StiC^RH  B  +  C<i>C(\R)Zh 
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Thus,  Equations  41  and  42  allow  us  to  convert  a  vector  in  one  coordinate 
system  to  another  coordinate  system. 

It  is  also  possible,  and  necessary,  to  express  the  angular  velocity 
(discussed  previously)  of  the  rotating  system  relative  io  the  fixed  system 
in  terms  of  the  rotating  triad  with  expressions  containing  the  Euler  angles  and 
their  derivatives.  Using  prior  notations  the  angular  velocity  can  be  written  as 


(ti 


dt  dt  dt 


(44) 


To  obtain  Z  in  the  rotating  coordinate  system,  i.e.  ZB,  determine 

-4  «*<+•♦ 

i"  in  terras  of  1B,  jg,  kg 

*4  4  •»  4 

j*  in  terms  of  iB,  jB,  kg 

■4  4  4  4 

kA  in  terms  of  iB,  jB,  k0 


Using  relations  already  established 


r  -  c^b-s^b 

k A  -  +  S*C<9jB  +  OfiCffkg  (45) 


Ii  like  components  are  now  combined  the  following  results  are  obtained: 


(w)xb  “ 

d<f> 

Us  e 

■4 

in  the  ig  direction 

dt 

dt 

(w)YB  - 

i£c* 

dt 

+  i^ScSC# 
dt 

in  the  j  0  direction 

(fjj)  - 

+  ^Lc<t>ce 

«* 

in  the  k  g  direction 

K  JZB 

dt 

dt 

(46) 
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It  is  also  possible  to  solve  for  and  jil  in  terms  of  (<.j)xb, 

dt  dt  dt 

and  (w)zb  obtaining 


(djj  g  +  myg  tan  6  S$  +  tan  ®  ^ 

6>yg  sec  0  +  uZBC^sec  9  (47) 

%e  now  have  sufficient  mathematical  tools  to  derive  the  equations  of 
motion. 

B.  The  Equations  of  Motion 

I.  Missile  Fixed  Coordinate  System 

In  this  lirst  derivation  of  the  equations  of  motion,  all  forces  and 
moments  are  referred  to  a  coordinate  system  rigidly  attached  to  the  mis¬ 
sile.  The  following  definitions  and  assignments  are  made  for  this  case: 


i£ 

dt 

d^ 

dt 

d 0 
dt 


ip  jj,  kj  A  space-fixed  (inertial)  coordinate  system  whose 
origin  is  located  at  the  center  of  the  earth  with 
the  kj  axis  coincident  with  the  spin  axis  of  the 
earth. 


* E’  i  e*  ^  coordinate  system  that  rjta^e^with  th$  earth, 

having  the  same  origin  as  ,  j  ,  k  j  with  k  £ 

-* 

coincident  withkj.  This  system  will  be  called 

the  rotating  earth  coordinate  system  and  is  that 
to  which  the  missile  trajectory  will  be  referred. 

iM,  jM,  kM  A  coordinate  system  rigidly  attache^  to  the  mis¬ 
sile,  origin  as  yet  unspecified,  but  i  taken  along 

the  longitudinal  axis  of  the  missile. 
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-*  4  4  -♦  H»  -4 

ftE  The  angular  velocity  of  (iE,  jE,  k£)  relative  to  (ip  jp  kj); 
that  is,  the  earth's  rotation. 

^  -4  4  -►  4  4  4 

a*  Angular  velocity  of  (iu,  jM,  kM)  relative  to  (ig,  jE,  kE). 

4  4  4,  4  4  4 

ft»F  Angular  velocity  of  (iM,  kM)  relative  to  (ip  jp  kj). 

®e  can  note  by  a  theorem  previously  proved  that 

-  -♦  + 

(Of  *  6)  +  Q  £  (48) 

-4 

Let  now  iu  be  along  the  longitudinal  axis  of  the  missile,  as  previously 

4 

stated.  Let  also  R  be  a  vector  defined  between  the  origins  of  the  inertial 

4  4  4  4  4  4 

coordinate  system  (ip  jp  kj)  and  the  missile  coordinate  system  (iM,  j  ^  k^). 

•4 

Finally,  let  SM  be  the  vector  from  the  origin  of  the  missile  coordinate  sys- 

4 

tem  to  the  gravity  center  of  the  missile,  and  R0  be  a  vector  from  the  origin 
of  the  inertial  coordinate  system  to  the  terminus  of  S||.  These  vectors  are 
shown  in  Figure  7. 


Figure  7 
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Ie  have  by  Newton's  Law  for  a  rigid  body 


d«*R 

IF  -m 


dt* 


(49) 


where  XF  is  the  summation  of  forces  acting  on  the  body  and 


dt* 
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the  acceleration  of  the  center  of  gravity  of  the  missile  relative  to  the 


inertial  coordinate  system.  In  particular,  the  expression  for 


dt* 


IS 


given  by  Equation  32  with  A  and  B  replaced  by  R  and  SM  respectively: 


dIh.  *£!“♦  2(aF)Mx  » ((3 ,»„» s„> .  (») 

dt*  dt*  dt*  dt  dt 


If  the  origin  of  the  missile  coordinate  system  is  located  at  the  missile 
gravity  center,  then  S  u  becomes  zero,  and 


IF 


d,*R0  d,*R 

a  — . .  •  t&  — -- 

at*  dt* 


(51) 


For  this  assignment,  the  fictitious  accelerations  (Coriolis,  etc)  in  Equa¬ 
tion  30  become  zero.  It  is  for  this  simplification  that  the  center  of  gravity 
of  the  missile  is  chosen  as  the  origin  of  the  "missile*  *  coordinate  system. 

To  obtain  a  particular  form  of  the  equations  of  motion  we  can  write 
d  R 

_i _  in  either  the  inertial  or  missile  coordinate  system  as  follows: 

dt 


d  R 

=  (V>1  *  (v)*i‘i  ♦  (V>Ylil  +  (V>z£l  =  Xlf  M  +  (VYM^M  +  ^VZM^M 

dt 


(52) 


Using  Equations  24  and  34,  we  can  express  Equation  51  in  the  follow' 
ing  forms 


IF  .  .  4Wl  .  .  400y 

<5  3a) 

dt  dt 

Ve  dy(V)  y  -• 

2.F  «•  m  +  (wf  )y  x 

dt 

<V)y 

(53b) 

•4 

*  VV>M 

EF  «■  m  - -  +  ((&>)  m  + 

dt 

*♦  -# 

(flE)M)  x  (V)y 

(53c) 

Note  that  Equation  53c  is  expressed  in  terms  of  the  missile  coordinate 
system.  To  use  this  equation,  the  forces  also  will  have  to  be  expressed  in 
terms  of  the  missile  coordinates. 


For  the  moment  equation  we  have 

v- .  Wh  ,  4<i>» 

dt  dt 


(54) 


there  IL  is  the  summation  of  torques  acting  on  the  missile,  and  (J^  -  (J)^ 

is  the  total  angular  momentum  of  the  rocket  relative  to  the  inertial  system, 
but  in  terms  of  the  ( y  fJt  or  (f^,  fy,  It^)  coordinates,  respectively. 

Again,  using  Equation  24 

*(*F)yX(j)y  (55) 

dt  dt 


which,  when  considering  the  earth's  rotation  again,  becomes 


y.(n  j„).l 

dt  dt  *4  t  m  H 


(56) 


If  the  principal  axes  of  the  missile  are  chosen  to  be  coordinate  axes  of  the 
missile  coordinate  system,  we  can  write  the  angular  momentum  as 

4 

!xu,  IYU,  Izu  being  the  moments  of  inertia  about  the  iy,  ju,  kM  axes, 

respectively.  Equations  53  and  56  are  the  general  vector  equations  of 
motion  of  a  missile. 


m 


M 
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Since  the  spin  rate  11 E  of  the  earth  is  known,  we  can  appropriately  se- 

*♦  -4  W 

lect  a  convenient  orientation  of  the  4,  jj ,  kj  and  iu,  jM>  kM  coordinate  sya- 

•4  «*  *• 

teas  and,  by  use  of  the  Euler  angles,  obtain  (1  E  in  terms  of  the  iM,  jM,  kM 
system,  that  is,  (n£)M. 

-a 

To  accomplish  this  objective,  set  both  kj  and  kE  permanently  along  the 
spin  asis  of  the  earth,  so  that 

*  **£*1  *  E  <5g> 


and  determine  the  components  of  k{  «  k£  in  terms  of  iM,  j  ^  k^.  Using 
Equation  42  yields 


«*  •* 
flgkg  «  — flgS&y  +  QgS^C^jy  +  flgCtfrCA  y  (59) 

In  component  form.  Equation  59  becomes 

<nE>  XM  * 

(fig)  Yk  *  ^^4>C8 

(fig  )  2||  ■  0  fC<j>C8  (60) 


— 

In  the  vector  Equations  53  and  56,  (unlikt  fig)  <if  and  V  are  not  known 

quantities.  They  are  unknowns  which  will  ultimately  give  rise  to  the 
trajectory  through  solution  of  the  differential  e  Rations  of  motion.  From 
before 


(V)M  -  Wxh'h  *  (V)Yu'u  +  (V)7uk 


zu  u 


(52) 


Similarly,  writing  Z  in  component  form  yields: 


«  (w)jfi^iy  +  (w)y^iM  +  IU 


(61) 
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Now,  we  can  immediately  write  the  basic  vector  equations  in  component 
form  which  are  valid  for  the  spherical  aad  rotating  earth  case.  Two  facts 
should  be  pointed  out  in  these  general  vector  equations:  (1)  Factors  such 
as  variable  air  density,  variable  wind,  and  the  gravitational  acceleration 
do  not  appear  in  the  inertial  terms  (which  have  just  been  defined),  but 
only  in  the  force  and  moment  part  of  the  equations  of  motion  (which  will 
be  derived  in  Section  8  of  the  PROCEDURE).  (2)  The  general  vector 
Equations  53  and  56  are  written  in  terms  of  the  missile  coordinate  system 

(as  indicated  by  the  presence  of  du(~)  ).  As  such,  the  unknown  qsanti- 

dt 

,  *  ■*  . 

ties  Vu,  cju  must  be  related  back  to  a  coordinate  system  to  which  the 
motion  is  referred.  The  latter  presents  no  additional  derivations.  RB  in 

Equation  43,  which  expresses  the  components  of  an  arbitrary  vector  hi  the 
coordinate  system,  may  be  replaced  by  the  missile  velocity  vector,  <V)M. 

Similarly,  (w)M  may  be  related  to  the  rates  of  change  of  the  Euler  angles 

by  Equation  46  or  47.  Vhen  the  angular  velocity  of  the  missile  is  known, 
the  new  Euler  angles  can  be  computed  by  an  integration,  while  knowing 
the  velocity  components  will  yield,  upon  integration,  the  missile  position 
( (R)jce*  <R)ye.  (Wze)*  These  equations  are  all  tabulated  at  the  end  of 
the  PROCEDURE.* 

i 

H.  Fixed  Plane  Coordinate  System 

It  is  not  necessary  to  specify  a  "missile  coordinate  system"  that 
is  rigidly  attached  to  the  missile,  as  has  been  done  previously.  In  this 
section  a  "missile  coordinate  system"  with  one  sxis  constrained  to  lie 
in  a  given  plane  is  considered,  although  a  coordinate  system  that  is 
rigidly  attached  to  the  missile  is  temporarily  used. 

The  following  definitions  and  assignments  sre  made: 


Same  coordinate  system  as  defined  earlier  (see  p  22) 
Same  coordinate  system  as  defined  earlier  (see  p  22) 


F«  the  interested  reader,  Reference  5  contains  the  derivation  of  these  equations  of 
notion  without  the  ose  of  matrix  notation. 
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*F‘  if*  k, 


*h» 


Qi 

>♦ 

n 


CJ 


The  missile  coordinate  system.  The  particular  orienta¬ 
tion  of  this  system  will  be  such  that  the  if  axis  is 

coincident  with  the  missile  axis  of  symmetry,  while 
the  jF  axis  is  constrained  to  lie  in  a  plane  parallel  to 

-*  +  **4  4 

the  plane  determined  by  ij,  jj.  The  origin  of  iF,  jF,  kF 
is  at  the  missile  center  of  gravity. 

A  coordinate  system  thst  is  rigidly  attached  to  the 
missile,  with  iM  along  the  missile  axis  of  symmetry, 

and  with  the  same  origin  as1F,  TF»^F* 

4  4  4  -4  *4  -a 

The  angular  velocity  of  ig,  j£,  kg  relative  to  ip  jj ,  kp 
Angular  velocity  of  ip,  jF,  kp  relative  to  iE,  jE,  kE. 

Angular  velocity  of  iM,  j^,  kM  relative  to  iF,  jF,  kf. 

Note  the  different  meaning  of  this  Z  compared  to  that 
defined  on  page  23. 


<0, 


ir  Angular  velocity  of  iM,  jM,  kM  relative  to  ip  jp  kj. 

Our  basic  vector  equations  in  terms  of  ip,  jF,  kF  now  become 
d*V, 


2F 


mJLL  +  «0)F  +  (Oe)*)  x  VF 

dt 


(62) 


1L  . 


dFJF 


+  ((0>F  +  (Oe^  *  Jf 


(63) 


In  these  equations  the  quantity  ((fi)F  +  (fig  )p)  is  the  angular  velocity  of 
the  missile  coordinate  system  to  which  we  refer  our  forces  and  moments, 

4  4  4  4  -4  4  ' 

relative  to  the  inertial  system,  i.e.,  (ip,  jF,  kp)  relative  to  (ip  jp  kj). 


To  derive  the  equations  of  motion  for^the^fixed  plane  system,  ic^is  first 
necessary  to  specify  the  orientation  of  if,  jF,  kF  relative  to  iE,  jE,  k£ 

in  terms  of  the  Euler  angles,  so  that  the  j  axis  lies  in  a  plane  parallel 


to  the  plane  determined  by  ij ,  jj.  This  condition  will  also  be  fulfilled  if 
the  jF  «uis  is  constrsined  to  lie  parallel  to  the  plane  determined  by  ig  •  ig< 

which  (by  definition)  is  coincident  with  the  ij,  jj  plane. 

To  accomplish  this  objective,  consider  the  following:  two  rotations: 

(1)  Rotate  theiE,  jE,  k£  coordinate  system  about  k£  by  an  angle  of  magni¬ 
tude  obtaining  the  identical  matrix  expressed  in  Equation  35.  (2)  Ro- 
(ate  the  resultant  coordinate  system  i',  j',  k'  about  the  j'  axis  by  an 
angle  of  magnitude  $,  the  final  position  being  denoted  as  i-.,  jp,  kp.  if 
the  resultant  matrices  are  multiplied  there  is  obtained 


-* 

*E 

C0Cd  -Stf  C^Sd 

*  F 

*E 

m 

sifrce  sv&sd 

-♦ 

-a 

kE. 

. -sd  o  cd 

“  kE  - 

(64) 


It  can  be  seen  from  the  inverse  that  jF  has  no  component  along  lc£  for  all 
\ft  and  d,  which  fulfills  the  constraint  on  the  jF  axis. 


The  angular  velocity  0  of  ip,  jF,  kF  relative  to  j£,  k£  can  be 
written  as 


q  .  £ikB  ♦  d±]' 


(65) 


dt  dt 

-♦  -*  ■+  “* 

By  obtaining  the  components  of  k£  and  j'  on  the  i£,  j  p  kF  axes,  we  have 


<fi> 


-set ±L  * 

dt  F  dt  F  dt  F 


(66) 


Solving  for  li!  and 
dt 


— 

dt 


yields 
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d  y 

m 

Wzt 

IT  * 

•in  9 

cca  6 

<o7) 


—  “  (fi>¥?  (68) 

dt 

For  Equations  62  and  63  to  be  of  use,  it  is  necessary  to  express  JF,  the 

angular  momentum,  in  terms  of  the  angular  velocities  and  the  moments  of 
inertia  of  the  missile,  namely 


Jp  -Ixf^F^Xf'f  +  ^F^F^Yf’f  +  IZF<<UF)2FkF  (69) 


where 

uip-w  +  n  +  Op  (70) 

Additional  relationships  between  and  fi  must  be  derived  since  six  un¬ 
knowns,  instead  of  the  usual  three  for  the  angular  velocity,  are  introduced 
_»  ■* 

when  m  and  fl  are  written  in  component  form.  These  additional  relation¬ 
ships  will  now  be  derived. 

It  will  be  recalled  that  both  iM  and  if  are  to  be  coincident  at  all  times. 
Hence,  the  only  rotation  enjoyed  by  the  iM,  jM,  kM  coordinates  not  shared 
by  the  ip ,  jF,  kp  system  is  the  motion  about  ip,  or  the  spin  of  the  missile 
about  its  axis  of  symmetry.  This  can  be  expressed  as 

(<u)XF  -  P  <">YF  -  (<o)ZF  -  0  (71) 

where  P  is  the  spin  of  the  missile. 

Since  a  new  unknown,  P,  has  been  introduced,  it  becomes  necessary  to 
eliminate  one  of  the  unknown  components  of  ?lF.  This  is  easily  done  since 

(fl)  can  be  written  in  terms  of  (fl)  by  use  of  Equation  67. 

X  F  £  r 


(fl)XF  =  “<n>ZF  ta°  9 


(72) 
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Finally,  we  cent  account  for  the  earth's  rotation,  flE,  by  defining  j,,  kj 

■nd  iE,  jE,  lcE  coordinate  systems  as  in  the  "missile  fixed"  case  with 
the  result  as  before 

*•  *♦*♦*♦- 
where,  by  Equation  64,  kE  can  be  expressed  in  terms  of  the  iF,  jF,  kj. 

system.  Therefore: 


(HeIxF  ” 

(73) 

(fiE>VF  “  0 

(74) 

(Oe)  ZF  *  (OE)C0 

(75) 

Kith  the  relations  obtained  from  the  above  analysis,  (»p)p  can  now  be 
written  as 

(cjp)p  ■  (P  -  tan  fl[fl)2p  -  (ftp  )S#)ip  +  (f^YE^F  +  W^ZF  *  (^E^ZF^^F  ^ 

Once  Equations  62  and  63  are  integrated  to  obtain  the  unknowns 
(v)xp  (V)YF,  (V)ZF,  p,  (0)yf»  (fl)zF’  the  inverse  th«  ®atr:s  in  Equa- 
tion  64  is  used  to  obtain  (V)XE,  (V)yE,  (V)zr.  Equations  67  and  68  are 

used  to  determine  JL t  and  i?.,  integration  of  which  provide  ^  and  8. 
dt  dt 

These  equations  are  tabulated  in  Part  V  of  this  section. 

HI.  Forces  and  Moments 

The  forces  and  moments  acting  on  a  missile  may  be  classified 
into  four  categories,  as  follows: 

1.  Gravitational 

2.  Aerodynamic 

3.  Jet 

4.  Guidance 
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Guidance  terms  are  beyond  the  scope  of  this  report  and  will  not  be  cod' 
sldered.  The  remaining  forces  and  moments  wilt  be  considered  in  the  or* 
der  shown  on  page  31. 

a.  Gravitational  Force 

The  gravitational  force  acts  a;  the  center  of  gravity  of  the  missile 
and,  hence,  does  not  produce  any  moments.  Further,  when  specifying  this 
term  and  other  altitude  dependent  terms,  distinction  must  be  made  between 
the  spherical  earth  and  the  flat  earth  cases  in  the  equations  of  motion.  In 
both  cases,  however,  the  magnitude  is  given  by  g  -  gjRw*/h0*  where 
is  the  mean  radius  of  the  earth,  and  g0  is  the  gravitational  acceleration 
at  sea  level.  To  ascertain  the  direction  of  this  force  for  the  flat  non-rotat- 
ing  earth  we  may  take  the  inertial  coordinate  system  ij,  jj,  kj  to  be  such 

that  k{  is  pointing  vertically  upwards. 

-*  •* 

The  gravitational  force  mgY  (Y  being  a  unit  vector  specifying  the 
direction)  then  becomes  mgY  -  -mgkj  (acts  opposite  to  kj).  With  respect 

•4  -*-+-» 

to  the  missile  fixed  system  (by  obtaining  kj  in  terms  of  iM,  jM,  kM),  this 
force  is 

mgY  «  —  S<pC#jy  —  C$C<?ky]  07) 

and,  by  similar  reasoning,  for  the  fixed  plane  system 

mgY  *  [S<?ip-C0kp]  (78) 

In  both  cases,  the  altitude  (hQ)  will  equal  (R)Z£. 

For  the  spherical  earth  case,  it  is  necessary  to  define  a  vector 
from  the  earth’s  origin  to  the  center  of  gravity  of  the  missile.  It  is  along 
this  vector  that  the  gravitational  force  acts. 
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Referring  to  Figure  8 


R-(R>Xe‘E  +(R)ye*E  +  Wzi^E 


R  -  RY 


Hence 


(79) 

(80) 


i  .1 

R 


(rWe 


WyeU 


v/(r)Je  ♦  <r4e  ♦  (*)•„  V(R)*xe  ♦  (*>»ys  ♦  WZB 


(R)zEkJ 


VW1  ”  <R)*  +  (R)1 

XE  YE  ZE 


(81) 


For  the  missile  fixed  coordinate  system, 

4-  CdC^iy  +  (S0S<£C^  -C<£S^)jM  +  (S8C^C^  +  S^S^)kM  (82) 


I 
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with  similar  expressions  for  jE  and  kE.  Combining  ail  the  iy  components 
will  yield  in  this  direction: 

(»fY)Xk  - - :2I - [(R)xeC5C0  +  (R)YEjCdS0  -<R)zeS03  (83) 

The  same  procedure  is  used  for  the  fixed  plane  coordinate  system,  where 
again  a  tabulation  is  presented  in  Part  V  of  this  section. 

The  altitude  for  the  spherical  earth  case  is  simply 


r  .  jRj  -  h0  -  V(R)*xe  +  <r>Ve  +  (r>ze  (SO 

b.  Aerodynamic  Forces  and  Moments 

The  aerodynamic  forces  and  moments  that  will  be  included  in  the 
equations  are  presented  ill  the  following  table.1  For  convenience,  both  the 
scalar  magnitude  and  direction  are  included  in  the  table. 


Syichal 

Fete* 

Mam ant 

Scalar  Magnitude 

Olractlan 

(DF) 

Axial 

Drag 

- 

P^DA 

Along  mis¬ 
sile  axis 

(NF) 

Normal 

(X^VrltN 

1  to  mis¬ 
sile  axis 
in  plane  of 

yaw 

(MF) 

Magnus 

- 

Pd*"FXVrkF 

1  to  plane 
of  yaw 

(RM> 

- 

Restoring 

pd'vr*k^(A-r) 

1  to  plane 
of  yaw 

(MM) 

— 

Magnus 

pd*eu  p  xVt  k  p(  A  p-r) 

_1_ 

i.  to  missile 
axis  in 
plane  of  yaw 

(YDM) 

Yaw  Damping 

pdW('c^Y  +<4z)2kM 

Independent 
of  yaw  plane, 
but  1  to  mis¬ 
sile  axis 

*  This  table  was  essentially  taken  from  Reference  6.  For  a  more  complete  listing  with 
explanations,  see  Reference  4. 
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Spnttal 

Fare*  Mawanf 

Scalar  Magaitarfa 

Otracftaa 

(AR) 

-  Roll  Damping 

',v'k*>  “tx* 

Along  mi  a* 

Moment 

aile  axis 

—  Roil  Moment 

pV,h^‘ 

Along  mis¬ 
sile  axis 

The  velocity  Vt  expressed  in  the  table  is  the  velocity  of  the  missile 

relative  to  the  wind.  For  example,  if  the  wind  has  a  velocity  component 
(VW)XM  in  the  same  direction  as  the  missile  velocity  component  (V)XM 

as  shown  below,  then  the  velocity  of  the  missile  relative  to  the  wind 
«Yt)xn)  becomes 


(V,)XM-(V)XM-(VW)XM 


(85) 


Simil  arl  y 

(Vf)yM-(V)yw-(Vt)yM 
lVz»l  “  (V)ZM  " 

It  should  also  be  noted  that  the  angular  velocities  «Fxt  ‘“py*  an<*  WF2 

expressed  in  the  table  are  the  three  components  of  *n  Tectof  Equation  48. 
For  example,  in  the  iy  direction,  with  the  use  of  Equation  59,  can  be 

written  as 


(<UF^XM  =  -(nE)S5  +  ^XM 


(87) 
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Finally,  if  the  wind  is  given  in  terms  of  either  the  inertial  or  rotating 
earth  coordinate  system,  the  Euler  transformations  should  be  used  to 
determine  the  wind  components  along  the  axes  of  the  missile  coordinate 
system. 

Each  of  the  forces  and  moments  given  in  the  table  will  now  be  briefly 
discussed.  The  reader  is  referred  to  other  texts  for  a  more  complete  dis¬ 
cussion. 

The  axis}  drag  acts  along  the  missile  longitudinal  axis  of  symmetry, 
but  in  direction  opposite  to  the  velocity  {-(Vf)Xy)  of  the  missile  along 

that  axis.  Since  the  line  of  action  of  this  force  passes  through  the  cen¬ 
troid  of  the  missile,  there  is  no  induced  moment  due  to  this  force.  The 
magnitude  of  this  force  is  as  given  in  the  table. 

The  normal  force,  on  the  other  hand,  acts  perpendicular  to  the  missile 
axis  and  lies  in  the  plane  of  yaw,  that  plane  determined  by  the  velocity 
vector  of  the  rocket(Vr)  and  the  rocket  axis  (See  Fig  U,  p  37).  The  com¬ 
ponents  of  the  normal  force  along  the  and  ky  axes  act  opposite  to  those 

of  the  missile  velocity  components  as  shown  in  Figure  10,  p  37.  To  ob¬ 
tain  the  magnitude  of  these  components,  as  shown  in  the  tabic  (pp  34  and 
35)  note  the  normal  force  lies  in  the  yaw  plane  but  perpendicular  to  the 
missile  axis;  hence,  if  y  (in  Fig  11)  is  defined  as  the  angle  between  the 
yaw  plane  and  the  jM  axis,  then  it  is  easily  seen  that 

|NF|ym  •  iNFl  cos  7 

|NpizM  -  |NF|  y  (88) 


where 


Cos  y 


<Vym 


sJ(\.Y  +(Vr)J 
V  ‘  YM  7.M 


(89) 


and 


Sis  y 


This  force  does  not  act  at  the  centroid  (CG)  of  the  missile,  but  at  what 
is  known  as  the  center  of  pressure  (CP),  which  is  determined  only  by  the 
exterior  configuration  of  the  missile.  Conse quendy,  this  force  gives  rise 
to  a  moment,  called  the  restoring  moment.  If  the  distance  between  the  mis* 
sale  center  of  gravity  and  the  center  of  pressure  (CP)  is  denoted  by  A'  (CP 
is  assumed  behind  the  CG  as  shown  in  Figure  12),  then  the  components 
of  the  restoring  moment  are: 

|RMjYM  -  |NFiZMA’ 

|RM!2W  -  !NFjYMA'  (90) 


Figure  12 

If  the  velocity  components  are  along  the  positive  j  andt:M  axes,  then, 

by  the  right-hand  rule,  the  components  of  the  restoring  moment  will  be 
negative  along  the  j  axis  and  positive  along  the  kM  axis. 

The  norma!  force  and  axial  drag  as  described  above  should  not  be  con¬ 
fused  with  the  lift  and  drag  which  are  components  perpendicular  and 
parallel,  respectively,  to  the  velocity  vector  of  the  rocket. 


The  magnus  force  depends  upon  the  angular  velocity  of  the  missile  sod 
sets  perpendicular  to  the  plane  of  the  yaw.  Consequently,  to  obtain  the 
component  of  this  force  along  the  jy  axis,  it  is  necessary  to  multiply  the 

scalar  magnitude  by  sin  y.  Thus, 

|MF|ym  -  |MF|  sin  y  along  jy 


and 


|MF|zm  “  |MF|  cos  y  along  iy  (91) 

If  (Vf) YMand  (Vt)  ZM  are  the  velocity  components  of  the  missile  relative 
to  the  wind  (discussed  previously),  then  the  sign  of  the  components  of 
this  force  is  as  indicated  in  Figure  13*  where  the  spin  is  considered  to  be 


<*>ZM 


Figure  13 

•m 

positive,  i.e.,  acts  in  the  direction  of  the  positive  iy  axis.  If  the  direc¬ 
tion  of  the  spin  is  reversed,  the  magnus  force  will  act  in  the  opposite 
direction.  Again,  since  the  magnus  force  does  not  act  at  the  centroid  of 
the  missile,  bin  at  the  magnus  center  of  pressure,  this  force  gives  rise 
to  a  moment  called  the  magnus  moment.  If  the  magnus  center  of  pressure 
is  again  behind  the  center  of  gravity,  as  indicated  in  Figure  14  (p  40), 
then,  by  the  right-hand  rule,  both  components  of  the  magnus  moment  will 
be  in  the  positive  jy,  kM  directions,  the  magnitude  of  which  will  be  the 

magnitude  of  the  scalar  component  multiplied  by  the  moment  arm  A  F  as 
seen  in  the  figure. 


39 


Figure  14 

The  yaw  damping  moment  tends  to  reduce  the  magnitude  of  the  yaw 
(assuming  kH  positive)  and  acts  opposite  to  the  resultant  cf  the  lateral 

angular  velocities  ((<uF)yM  (wp)ZM)-  This  is  shown  in  Figure  15.  Note 

that  there  is  no  component  of  the  yaw  damping  moment  along  the  missile 
axis  of  symmetry.  Noteworthy  also  is  that  a  force  associated  with  kH 

exists;  because  it  is  difficult  to  measure  this  effect,  however,  this  force 
was  omitted  in  the  analysis. 


(YDM) 


Figure  15 

The  final  aerodynamic  terms  to  be  considered  are  the  roll  damping  mo¬ 
ment  and  the  roll  moment.  Both  moments  act  along  the  longitudinal  axis  of 
the  missile.  The  former  is  induced  by  the  fin  cant,  while  the  latter  tends 
to  reduce  the  axial  spin  by  skin  friction  on  the  surface  of  the  missile  dur¬ 
ing  flight.  Both  moments  act  in  a  direction  opposite  to  that  of  tht  spin  of 
the  missile. 
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It  should  be  noted  that  none  of  the  aerodynamic  coefficients  are  assumed 
to  be  linear  in  nature  and  hence  are  not  to  be  taken  as  slopes  from  experi¬ 
mental  curves.  Rather,  these  coefficients  are  obtained  directly  from  experi¬ 
ment  as  functions  of  Mach  number  and  angle  of  attack.  For  example,  one 
value  of  Mach  number,  the  angle  of  attack,  and  the  corresponding  axial 
drag  are  sufficient  *o  determine  one  value  of  kcA. 

c.  Jet  Forces  and  Moments 

The  jet  forces  and  moments  considered  in  this  report  are  due  to 
thrust,  thrust  malalignment  torques,  and  jet  torque  to  initially  spin  the 
missile  about  its  axis  for  spin  stabilization. 

The  jet  thrust  is  the  primary  force  that  imparts  forward  motion  to  the 
missile.  Although  it  should  ideally  act  along  the  missile  axis  of  symmetry, 
the  thrust  may  have  components  along  all  three  axes  of  the  missile  co¬ 
ordinate  system.  This  is  due  to  imperfections  in  the  rocket  design,  the 
practical  nature  of  propellants,  and  other  factors.  To  obtain  quantitative 
relationships,  one  may  conveniently  define  two  angles:  (1)  ST,  the  angle 

between  the  direction  of  the  thrust  and  an  axis  parallel  to  the  missile  axis, 
and  (2)  SA,  the  angle  between  the  projection  of  the  thrust  force  on  the 

jM,  kM  plane  and  the  jM  axis.  These  angles  are  illustrated  in  Figure  16. 


Figure  16 

From  geometrical  considerations. 


(T)xm  =  T  cos  8r 

along  iM  axis 

<T>VU  -T..S,  SA 

-* 

along  jM  axis 

(T)zfl  »  T  sin  5-j-  sin  S  A 

-* 

along  kM  axis 

(92) 
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If  the  direction  of  the  thrust  does  not  change  relative  to  the  missile, 
then  for  the  missile  fixed  coordinate  system,  both  &A  and  $T  can  be  taken 

as  constants  during  the  missile  flight. 

To  obtain  the  corresponding  malalignment  torques,  reference  is  made 
to  Figure  17. 


If  (T)XM,  (T>ym.  (T)zm  act  alon8  the  positive  iM ,  jM,  kM  axes  respectively, 
then  taking  moments  will  yield 

Thrust  moment  in  iM  direction: 

(TM)xm  -  (T)zM*va  “  (T)YMrl  (93) 

Thrust  moment  in  jM  direction: 

(™>YM  =  <TWl  "  <TWt  <94) 

Tlirust  moment  in  direction: 

(TM>zm  =  (T)YM't  -  (T)XMru  (95) 

where  r  ,  r. ,  and  r  are  defined  in  Figure  17. 
i  ’  1  u 
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Finally,  the  jet  torque  can  be  considered  as  proportional  to  the  total 
thrust  or  KT.  If  the  spin  about  the  i^  axis  is  clockwise,  then  the  jet 
torque  is  said  to  be  positive  by  the  right  hand  rule.  Its  purpose  is  to  ini* 
daily  spin  the  missiie  for  spin  stabilization.  It  acts  as  indicated  in  Fig* 
ure  18. 


Figure  18 

In  summary,  the  forces  and  moments  as  discussed  in  this  section  are 
exhibited  in  the  table  shown  on  page  44. 

IV.  Initial  Conditions,  Conversions  to  Specific  Cases  and  Singulari* 
ties 

a.  Initial  Conditions 

Initial  conditions  for  the  equations  of  motion  will  consist  of 
initial  values  for  the  following  variables:  (R)x  (R)YJ  .  (R)Zp  (V)XM, 
(V)ym.  (V)zM;  0,  <Jj,  <£;  (u)xMi  (cu)ym,  («)zm  •  Of  these  quantities  (R)xi, 
(R)yi«  (R)Z|,  and  6,  «4,  «j>  are  in  reference  to  the  inertial  coordinate  system, 
while  (V)XM,  (V)VM,  (V)ZM  and  (cj)xm,  ^ zh  are  rffferrcd  t0 

missile  coordinate  system.  In  addition,  it  is  important  to  know  the  orienta¬ 
tion  of  these  missile  coordinates  at  launch,  relative  to  the  earth.  Any  ini¬ 
tially  directed  forces  (for  example,  thrust  malalignment)  will  influence  the 
particular  latitude  and  longitude  of  the  missile  impact  point  at  the  end  of 
the  trajectory.  These  directed  forces  take  on  added  importance  in  influ¬ 
encing  the  terminal  point  of  the  trajectory  if  an  ellipsoid  is  used  to  ap¬ 
proximate  the  earth’s  geometry. 
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Components 


•rust  Malalignment  T  sin  ST  sin  SAru  T  cos  T  sin  STcos  SAt, 

Torque  -T  sin  5T  cos  SArj  -T  sin  SA  sin  $j*f  -T  cos  5Ar„ 


■  <a«. 


To  obtain  the  initial  conditions,  let  us  now  precisely  orient  the  inertial 
coordinate  system.  As  previously  stated,  the  origin  of  this  coordinate  sys- 
tern  is  situated  at  the  earth’s  center  with  the  kj  axis  coincident  with  the 

a*is  of  spin  of  the  earth.  Now  let  ij  (arbitrarily)  coincide  with  the  Green- 
wich  Meridian  (zero  degrees  longitude),  while  jj  is  perpendicular  to  ij 
and  kj  in  a  direction  given  by  the  right-hand  rule.  Further,  let  us  assume 
the  existence  of  another  coordinate  system  L.,  jc,  kc  initially  coincident 
with  ij,  jj,  kj.  Upon  successive  specified  rotations,  it  is  intended  that 

*C»  ic*  become  the  missile  coordinate  system.  Here  the  subscript 

C  simply  denotes  that  this  coordinate  system  is  used  to  determine  initial 
conditions.  As  before,  primes  will  be  used  to  denote  the  intermediate 
rotations.  Let  the  launch  point  of  the  missile  be  at  A°  longitude  and  B° 
latitude. 

Out  first  objective  is  to  orient  the  ic  axis  so  that  it  passes  through  the 

launch  point  of  the  missile.  This  is  easily  done  by  two  rotations,  first  a 
rotation  about  k(  of  magnitude  A°,  and  second  a  rotation  about  jc*  of 

-*  «♦  4  4 

magnitude  B°.  Note  that  i^;,  jc,  k  -  is  now  labelled  as  ic''*  ic"»  kc"  • 
Using  the  techniques  in  Part^A  of  the  Procedure,  we  can  readily  establish 
relationships  between  ij,  jj,  kj  and  »c",  jc",  ^c"'  ®y  matrix  multiplica* 

tion,  we  have 


- 

-♦ 

1  < 
*c 

CACB  CBSA  -SB 

-* 

ic” 

«■ 

-SA  CA  0 

h 

-* 

•* 

kc"  j 

SBC  A  SBSA  CB  _ 

.ki  - 

(96) 


It  will  be  convenient  if,  in  principle,  we  now  translate  the  i£ ’, 

coordinate  system  out  to  the  surface  oi  the  earth  so  that^the  origin  is  at 
the  launch  point  O'  of  the  missile.  It  will  be  noted  that  i^"  is  still  in  a 

direction  extending  radially  outward  from  the  center  of  the  earth  and  pro¬ 
vides  a  perpendicular  line  to  the  earth’s  surface  at  the  launch  point, 
while  j*c"  and  £c"  determine  a  tangent  plane  to  the  earth's  surface  at  O'. 
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Note  also  that  the  k  c  *  *  axis  points  towards  north.  In  this  "translated” 

position  O',  ic”»  jc"»  kc"  will  be  called  the  geographic  system.  To 
completely  specify  the  orientation  of  the  missile,  two  additional  angles 
are  specified  relative  to  this  geographic  system.  In  this  analysis  we  will 
use  (1)  the  angle  of  declination  of  H6  from  the  perpendicular  to  the  earth's 
surface,  and  (2)  a  lateral  orientation  angle  G°  measured  from  the  nejptive 
kc"  axis  in  a  counterclockwise  direction  in  the  tangent  plane. 

These  angles  are  pictured  in  Figure  19. 


Narth 


Figure  19 
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The  object  now  is  to  perform  two  additional  rotations  using  the  input 
information  of  G°  and  H°,  so  that  the  i~"  axis  finally  coincides  with  the 

missile  axis.  We  first  rotate  about  ic"  by  G°  so  that  -k  c"  coincides  with 

the  projection  of  the  missile  axis  on  the  jc",  k  c"  plane-  The  second 

rotation  is  about  the  axis  of  magnitude  H°,  which  finally  causes  the 

ic,M  axis  to  become  coincident  with  the  missile  axis.  It  is  important  to 

note  that  the  jM  axis  is  located  in  the  tangent  G°  from  the  east,  measured 

in  a  counterclockwise  direction.  Figure  20  shows  the  final  rotation. 


Mathematically,  the  resultant  matrix  for  the  G*  and  H  rotations  is 


lU 

CH 

SH.'IG 

-SHCG 

-* 

•  Jl 

*c 

■4 

m 

0 

CG 

SG 

>c' 

kii 

.  SH 

-CHSG 

CHSG 

-♦ 

*c 

(97) 


The  matrix  for  all  four  rotations  becomes,  by  matrix  multiplication. 


•* 

r  • 

'  -  1 

lu 

A„ 

*1 

-+ 

h 

-b 

kM- 

. 

lki  J 

where 

Au  -  CACBCH  -  SASGSH  -  CASBCGSH 

AtJ  -  SACBCH  +  CASGSH  -  SASBCGSH 

A„  -  -CHSB  -  CGCB5H 
A,  i  -  -SACG  +  SGSBCA 

(98) 

Alt  •  CGCA  +  SGSBSA 
A,  |  “  SGCB 

A„  -  SHCBCA  +  CHSGSA  +  CHCGSBCA 

Am  -  SHCBSA  *  CHSGCA  +  CHCGSBSA 

A„  «  -SHSB  +  CHCGCB 

We  are  now  in  a  position  to  obtain  the  initial  Euler  angles,  so  that  i{ 
becomes  coincident  with  the  missile  axis.  We  first  perform  a  rotation  about 

-4  *• 

kj  of  magnitude  ift,  and  a  rotation  about  j*  of  magnitude  6.  We  will  now  have 
-0  ■* 

i'"  coincident  in  both  cases.  To  obtain  jM*  we  perform  a  rotation  about 
i'"  of  magnitude  Mathematically,  the  three  rotations  expressed  in  matrix 
form  are 


48 


•u 

cdC4> 

cesy, 

-se 

-* 

m 

-  c4>s>p 

S&S$Sip  +  CipCtfa 

s^ce 

«* 

_  &8C<£>CiJj  +  S4>s^ 

$ec<f>s -  s^»c^ 

c<K& „ 

(99) 


This  and  the  preceding  matrix  must  be  equivalent,  since  the  final  positions 
of  the  missile  coordinate  systems  are  to  be  identical.  To  obtain  6,  <p,  we 
simply  compare  corresponding  elements  of  the  two  matrices.  In  particular, 


S6  -  SBCH  +  CBCGSH 

»\i*Cd  .  (SACBCH  +  CASGSH  -  SASBCGSH) 

*  ran  y>  *•  ' _ _ 

Q4>ce  (CACBCH  -  SASGSH  -  CASBCGSH) 


.  tan  <t>  *  SGCB- 

cec<t>  CHCGCB  -  SHSB 


(100) 


It  is  to  be  understood  that  these  angles  are  at  t  -  0,  or  at  time  of  launch. 

One  must  be  careful,  however,  in  ascertaining  the  correct  quandrant  of 
these  angles. 


For  the  initial  position  of  the  missile,  only  the  latitude  and  longitude 
are  required  to  compute  the  (R)XJ,  (R)YI,  (R)ZI  coordinates  of  the  missile 

launch  point,  fe  can  assume  at  t  «  0  that  the  reference  coordinate  sys¬ 
tem  rE,fE,tEis  coincident  with  the  ij,  jj,  lj  coordinates,  where  again 

the  ij  axis  coincides  with  0°  longitude.  From  Figure  21  (p  50)  it  is  easily 
seen  that 


(R)XI  “  %  co.  B°  co.  A0 

(R)yi  -  R»co.B°ainA°  (101) 

(R)z|  -  Rv  ton  B° 

Rw  being  the  mean  radius  of  the  earth. 

Similarly,  knowing  the  coordinates  of  the  terminal  position  of  the  mis¬ 
sile  enables  the  determination  (by  the  same  equation)  of  the  latitude  and 
longitude  of  the  trajectory' s  end. 
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Figure  21 

Finally,  knowing  the  orientation  of  the  jM  axis  and  any  directional 

forces  and  moments  relative  to  a  specified  compass  direction  enables  one 
to  compute  the  components  of  these  forces  and  moments  along  the  missile 
coordinate  system.  For  example,  if  one  wanted  the  effect  on  the  terminal 
point  of  the  trajectory  due  to  thrust  malalignment  initially  directed  in  an 
eascerly  direction,  then  d A  (the  lateral  thrust  malalignment  angle)  may  now 

be  given  a  definite  value,  namely  the  G°  as  we  have  defined  it.  Similarly, 
knowing  the  orientation  of  the  missile  coordinates  will  enable  the  compu* 
ration  of  any  initial  lateral  velocities  and  angular  velocities.  This  com¬ 
pletes  the  discussion  on  initial  conditions. 

b.  Conversion  to  Specific  Cases 

For  short  trajectories,  generalization  to  the  "rotating  earth"  and 
"spherical  earth"  cases  may  not  be  warranted  in  solving  the  equations  of 
motion.  Converting  the  equations  to  the  flat  earth  requires  using  the  appro¬ 
priate  gravitational  term  as  well  as  setting  the  altitude  equal  to  Z  (in 

place  of  ho  -  VWJce  +  ^ye  +  (r)*ze^  To  convert  t0  che  non-rotating 
earth  it  is  only  necessary  to  set  fiE  equal  to  zero.  Similarly,  to  neglect 
wind  one  sets  Vw  equal  to  zero. 

c.  Singularities 

Care  must  be  taken  that  certain  quantities  do  not  become  increas¬ 
ingly  large,  without  bound,  in  solving  the  equations  of  motion.  Such  a 
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singularity  is  present  in  the  missile  fixed  system.  Equation  47  is  repeated 
below  for  convenience 


(106) 


d^>  da 

-  c  __  —ft) 

dt  dt  X 

which  is  also  obtained  from  Equation  46  with  $  »  90°.  If  in  addition  to 
8  »  90°,  both  <uy  and  <az  equal  zero,  recourse  must  be  taken  to  higher  or* 

der  derivatives. 

V.  Tabulated  Equations  of  Motion 

Following  is  a  summary  of  the  si  x-degree-of -freedom  equations  of 
motion  for  body  fixed  system  (rotating  spherical  earth). 

Fore*  Equations 


T  co*  Sx-pd’( Vt)MkDA  +■ 


v'<R>'xE*(B)'vE*'R)iE 


t-CK)XE  coo  8  cos  0  -  (R)y£  cos  8  • in  0  +  (R^g  “b  ^  “ 
d(V)vu 

ra  - ((<u)VM  +  (flE)  sin  #  cos  8)  (V)ZM 

dt 

”  +  Cos  4*  co*  OMyii] 


T  sin  ST  cos  5a  - 


P^VmM^YM 


i/<v,)-y||  +  0*>*Zi( 


pd*(<uF)xM(Vt)MlcF(Vr)ZM 
\f(-\ )yM  +  (Vf  ) ZU 


V(R>X  E  T  ^YE  +  (R^ZE 


[-<R)  x^sin  $  sin  4>  co*  0  -  cos  <£  sin  0) 

-  (R)yE(sio  #  sin  0  sin  0  +  cos  0  cos  <£)  -  (R)ZE(sin  cos  0]  = 
f  d(V)VM 

m  I  - +  ((^ZM  +  (^E^  cos  ^  cos  ^  (^)xm 

l  dt 

-  «<u)XM  -  (fi£)  sin  $)  (V)ZM 


'1? 


T  sin  5t  sin  SA  -  pd\\  >kM\  >ZII 


^  >yh  ♦  <v,  >z»i 


pci'Cw F ) X Vt  ) F(Vr )  Y m 


mg 


v<v  )'vu  *  <V  >‘zu  vWXE  *  <R>Ve  *  <«>i£ 

[-(R)xe  (sin  6  cos  tf>  con  +  sin  <+>  sin  ) 

-(R)  y  E  (sin  6  cos  $  sin  0  -  sin  ^  cos  —  (R)zE(CoS  4>  cos  #)  ] 


m 


d(V)ZM 

dt 


+  ((^)XM-(nE)  sin  «(V)YM 

-  «^)VM  +  (fiE)  «n  4>  cos  0)  (V)XM 


Momont  Equations 

T  sin  5t  sin  5Atu  -  T  sin  5j  cos  5Arj  +  kT 
*  P(Vf)|Jc.£p(wF)xMd*  -  p(Vr)Mk?id*  - 


1  XM 


- -  (ftg)  cos  6  cos  <i>  (w)yy  *  (fi£)  cos  ^  0  (<«*)  ZM 

dt  ‘  '• 

+  f  1  ZM‘  1  YM  1  ((W)YM(CU)ZM  +  <nE>  si°  *  C°*  *(w)ZM 
+  (Q£)  cos  4>  cos  #(w)ykj  +  (0E)*  cos*  5  sia  <f>  cos  <£) 


T  cos  8-j.fj  -  T  sin  8  A  sin  STrt  - 


^Wn(A-,)(v.>zw 


v(V)»  +  (V  )J 

'  r  YM  t  ZM 


pd\ajF;XM(  V  )Mkp(  Af  -  r)(v  )  VM 


-  kH(V)pd^F)YM  - 


^'^YM  +  ^V,)ZM 


I  ym  |  +  (fig)  cos  0  cos  $  {«>XM  +  (^e)  *i°  $  (w)7M  j 

L  dt  J 

+  fl  XM~  I  ZkJ  (MZ«{«)XU  ~  (^E)  ®in  #  <»  ZM 
+  (fig)  cos  4>  cos  6  wx  -  (0 co»  0  cos  0  sin  9) 


T  sin  ST  cos  SAct  -  T  cos  5yru  + 


^'(V  )^kNa-f)  (Vf)VM 


V<vf>Vu  +  (vt)*ZM 


+ ^  ( v,  >**<*  F)  ZM 

V<Vym  +  <X)‘zM 

d(cu) 


lZii 


1 ZM 


dt 


-  (fiE)  sin  0  cos  G  (w)xk.,  -  (^e>  s‘“  9  (cj)  yM 

+  ^  YM“  !xJ  ^^Xi/^YM  +  *in  ^  co*  XU 

-  (fig)  sin  Q  (cj)ym  _  (Qe)‘  *in  9  sin  0  cos  9) 


Additional  Equations 

VXE  "  co*  ^  co*  ^  (V^XM  ■*■  (*°  ®  0  cos  ^  “  co*  0  *‘a  (v)  yM 

+  (sin  G  cos  0  cos  0  +  sin  <p  sin  0)  (V)ZM 

Vyg  «  cos  6  sin  0  (V)x«  +  (*in  6  sin  0  sin  0  +  cos  0  cos  0)  (VJyy 

+  (sin  6  cos  0  sin  0  ~  sin  0  cos  0)  (V)ZM 

VZ£  -  -sin  6  (V)  xm  +  sin  0  cos  <3  (V)yM  +  cos  0  cos  G  (V)  Zm 


*  ^VI^XE  co*  ^  co*  ^  +  ^VW^YE  co s  $  sin  0  —  (VW)Z£  sin  0 

(V^)  YU  “  (vv)  xE  **‘n  ^  ,in  0  co*  0  ~  cos  0  sin  0)  ♦  (Vw)yE 

(sin  0  sin  0  sin  0  +  cos  0  cos  0)  +  (V^ZE  *‘n  ^  cos  ^ 

('Vzw  "  *V»^XE  ®  co*  0  cos  0  +  sio  0  sin  0  )  +  (Vw)yE 

(sin  5  cos  0  sin  0  -  sin  0  cos  0)  +  (v^)Ze  Co*  ^  co*  ^ 
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<VA*  “ 

(Vym  " 

(Vf)ZM  . 

Equation* 

•in  G(0) 
too  0(0) 

too  0(0) 
coo  8(0) 

(R)XI(0) 

(R)yi<0) 

(R)ZI(0) 


d$ 

dT 


te 


OJy  COB  0  -  <0Z  *in  0 


d0 

dt 


UX  *  UV  t*’1  ^  **n  0  +  ^2  tafl  ^  coa  0 


d0  , 

_  «  aiy  sin  0  sec 

8  +  <U2  coo  0  sec  $ 

dt 

<V>XI|-<Vxu 

(V>yM  ~  (VV)  YM 

(tyF)YM  “  (a>)y m  +  **“  ^  co*  e 

(V)ZM-(VW)Zii 

^F^ZM  “  <«) Zk  +  co*  0  co*  ^ 

for  Initial  Condition* 

m  cob  H  sin  B  +  *in  H  co*  G  co*  fi 

_ sin  G  coo  B _ 

co*  H  coo  G  cos  B  -  oin  H  oio  B 

-  cos  H  cos  B  ota  A  +  sin  H  sin  G  cos  A  -  oin  H  cos  G  oin  B  oio  A 

_  coo  H  co*  B  sin  A  +  oio  H  oio  G  coo  A  -  oin  H  co*  G  oio  B  oio  A 

oin  0(0) 


-  Rv  cos  B  cos  A 

_  ...  if  liven  longitude  sad  latitude 

•  Rv  cos  B  mn  A 

-  Rv  sin  B 
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Stop  Equation  * 


(r)xe  = 

A 

I 

<V>XE* 

e  = 

f 

iidt 

0 

O 

dt 

(R)yE  - 

c 

(V)YEdt 

<6  = 

S 

i 

'o 

i£dt 

i it 

(R)  ZH  = 

c 

(V)zpdt 

t 

fo 

iiidt 

dt 

Equation*  for  Singular  Condition* 

For  0  -  90°  (<u)yM  /  0 

(&)  ZM 
dt 


d(ai)  yM 


dt 


)  ”  “>VM  l“—  I 


(*)* 


YM 


cosJ  <t> 


For  0  =  9Qa  (oj)ym  =  0  («)ZM  ^  0 


d<£ 

dt 


(w) 


ZM 


dt 


'  dt 


(oj)1 


ZM 


sin*<£ 


In  both  cases 


dy  _  d<6  ,  . 

-  *  -  -  (t u )  XM 

dt  dt 


Following  is  an  alternate  set  of  more  important  equations  for  the  fixed 
plane  coordinate  system: 


IF 


XF 


m 


dV 


XF 

dt 


IFyf  ^  m 


dV 


YF 


dt 


*nYFVZF-(!!ZF*inE,Cl)’VVF 


*  m  ZF  *  V®  VXF  -  <-n  ZF  •*"  °  -  "e*”  VZF 


S6 


ifzf  -  » 


dV 


ZF 


1  di 


+  ( ~flZF  ~  S5)Vyp  ~  (  0  y  pJVjjp 


^LXF  -  *XF 


dp 


~ft£C0S0flyp  -  i£il£  tta  0  -  Azf^yf  »«c*  <? 

dt  dt 


fiyF(  ^ZF  "*■  ^^ZF  “  (  f^zp  +  5) Oypiyp 


dO 


2LyF  -  - 1£-I  YF+  (nZF  +  nECe)(P-flzF*«  ^  “  GfiS#)!  XF 


dt 


-  (-0  Zp  t»  $  -  QgS  $)  ( nZF  +  flg  C  0)IZF 


2L 


ZF 


dn 


ZF 


dt 


*•  Q^S  50yp 


IZF  +  (  ”^ZF  **°  ^  “  fig®  ¥F 


“(fiyp)  (P  “0ZF  tWJ  tf-OgS*)!^ 


VXE  *  CltC$VXF  -SviVyp  +  C^S0VZF 
VyE  -  S^C«VXF  +  QfrVyjf  +S^S«VZF 

vZE  *  -sevXF  +  cevZF 


(VW)XF  -  (Vw)XEC0C$  +  <Vw)YES^C0-(Vw)ZES* 
<V*)yfr  -  -<V,)XES^+<VW)yEC^ 

(VW)ZF  -  (V1)yEC0S0+<Vt)yES^S0*(V1)ZEC« 


il  -  flyp 

dt 


— -  *  ^ZF  *ec  0 
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Flat  Earth 


-t  -*  -* 

mgy  =  reg  {  S0i  f  -  Ct/kj,  3 

Spherical  Earth 

mgy  -  — - B1 . .  ——  [  (-(R)xeC0C@  -  (R)yeS^C@ 

vCR)xe-<r4e  *<R>ze 

+  (R)ZES3)rF  +((R)XES0-(R>y.EC^)rF 

-* 

*  ((RiXEC0Sd-(R)YEStfS0  +  (R)ZE)kp  ] 

RESULTS  AND  DISCUSSION 

The  complete  set  of  the  equations  of  motion  are  tabulated  in  the  pre¬ 
ceding  section.  Before  these  equations  are  solved,  however,  the  user  must 
specify  orher  condtrions  of  the  problem.  For  example,  the  time  to  the  cut¬ 
off  point  (i.e.,  where  the  thrust  period  terminates)  must  be  given.  Varia¬ 
tions  in  mass,  cencer  of  gravity,  center  of  pressure,  and  magnus  center  of 
pressure  during  and  after  the  thrust  stage  are  required.  The  thrust  mal¬ 
alignment  distances  (rt  ,  rj ,  ru )  of  the  rocket  must  also  be  furnished.  Fur¬ 
ther,  a  complete  set  of  the  aerodynamic  coefficients  must  be  given  for  the 
complete  velocity  range  and  for  all  angles  of  attack  encountered. 

If  it  is  desired  to  include  wind,  a  wind  profile  must  be  available.  In  addi¬ 
tion,  air  density,  speed  of  sound,  and  gravitational  acceleration  as  func¬ 
tions  of  altitude  must  be  specified.  The  user  should  also  furnish  informa¬ 
tion  to  determine  the  end  of  the  trajectory  such  as  terminal  altitude,  etc. 

Similarly,  the  user  must  decide  upon  a  particular  numerical  integration 
scheme  in  order  to  solve  the  equations  of  motion.  Frequently  used  in  this 
area  is  one  of  the  Runge-Kutta  techniques. 

It  should  also  be  noted  that  in  the  tabulated  equations  of  motion  (for  the 
missile  fixer*  system)  several  terms  contain  a  trigonometric  function  of  the 
angle  <p.  For  spin-scabiliaed  rockets,  i.e.,  large  (w)  and  (<oF)XM, 

the  rate  of  change  of  may  be  quite  high  (d<$/dt).  This  may  require  small 
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cimc  increments  far  the  numerical  solution  of  these  equations  to  obtain 
suitable  accuracy.  This  problem  is  not  encountered  when  one  considers  a 
missile  triad  that  does  not  spin  with  the  missile.  For  this  latter  case, 
however,  one  seems  to  be  confronted  with  time  variant  moments  and 
products  of  inertia  for  asymmetric  missiles. 

A  final  word  about  the  equations  is  that  there  is  no  estimate  of  the 
dispersion  of  the  missile.  This  requires  the  computation  of  several  tra¬ 
jectories,  each  for  slightly  different  initial  conditions,  and  appropriate 
statistical  combinations  of  the  various  range  deflections  from  the  stand¬ 
ard  range. 
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APPENDIX  -  REQUIRED  INPUT  PARAMETERS 


The  complexity  of  the  six-degrce-of-freedcm  equations  of  motion  requires 
the  use  of  a  computer  for  solution.  These  equations  have,  therefore,  been 
programmed  for  the  IBM  709  of  Picatinny  Arsenal.  This  appendix  tabulates 
the  input  data  tnat  the  user  must  supply  to  the.  computing  personnel  for  a 
meaningful  statement  of  the  problem. 

Tab!*  of  Input  Dots 


1.  Specify  flat  or  spherical  earth  case. 

2.  Specify  rotating  or  non-rotating  earth  case. 

3.  Specify  phases  of  flight  that  are  to  be  considered. 

Phase  I,  Acceleration  of  booster  and  main  stage. 


Phase  II.  Coasting  of  booster  and  main  stage. 
Phase  111.  Separation  of  booster  and  main  stage. 
Phase  IV  Coasting  of  main  stage. 


Phase  V.  Acceleration  of  main  stage. 

Phase  VI.  Free  flight  of  main  stage. 

4.  Specify  the  time  duration  of  each  phase  considered. 

5.  Specify  mass  of  missile  and  booster  combination. 

6.  Specify  mass  of  booster  alone  (at  launch). 

7.  Specify  rate  of  change  of  mass  of  booster  during  Phase  I. 

8.  Specify  rate  of  change  of  mass  of  main  stage  during 
Phase  HI  (if  considered). 

9.  Specify  rate  of  change  of  mass  of  main  stage  during 
Phase  V  (if  considered). 
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iU.  Give  data  relating  thrusi  vs  time  for  Phase  I  and  for 
Phases  ill  and  V  (when  considered). 

11.  Specify  initial  moments  of  inertia  of  booster  and  main  stage 
(at  launch.  Phase  i)  along  principal  axes  of  rocket, 

12.  Specify  rate  of  change  of  moments  of  inertia  oi  booster  along 
principal  axes  during  Phase  1. 

13.  Specify  rate  of  change  of  moments  of  inertia  of  main  stage 
during  Phases  III  and  V  (when  considered), 

14.  Specify  distance  from  missile  nose  to  missile  centroid  at 
start  of  Phase  1  and  after  separation  of  booster. 

15.  Specify  rate  of  change  of  missile  centroid  position  for 
Phase  [  and  for  Phases  31!  and  V  (when  considered). 

16.  Specify  same  information  as  in  Step  14  for  missile  center 
of  pressure. 

17.  Specify  magnus  center  of  pressure. 

18.  Specify  maximum  diameter  of  missile  for  Phases  1  and  III. 

19.  Specify  magnitudes  for  the  aerodynamic  coefficients 

kDA  kd> 


k  F  kH 

as  functions  of  Mach  number  and  angle  of  attack  for  booster 
and  main  stage  combination  and  main  stage  alone. 

20.  Specify  thrust  malalignment  angles  (<5X,  §A)  for  Phases  I,  111, 
and  V. 

21.  Specify  thrusi  malalignment  distances  and  their  rates  of  change 
for  Phases  I,  III,  and  V. 

22.  Specify  proportionality  factor  of  thrust  used  for  inicir'  spin 
stabilization. 
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23.  Specify  latitude  and  longitude  of  launch  point  of  missile. 

24.  Specify  orientation  angles  of  missile  relative  to  geographic 
system. 

25.  Specify  initial  velocity  components  of  missile  ((V)XM,  (V)YM, 
(V)ZM). 

26.  Specify  initial  angular  velocity  components  of  missile 

27.  Specify  wind  velocity  components  as  function  of  altitude. 

28.  Specify  altitude  at  which  trajectory  is  co  terminate  (air  burst, 
ground  burst). 

29.  Specify  output  data  desired  ({R)x  £,  (R)yE  .  (R)ZE  >  0.  <£> 
(V)XE...). 

30.  Specify  time  increments  for  each  phase  when  above  quanti* 
ties  are  to  be  tabulated. 

31.  Specify  any  peculiarities  of  the  missile  system. 

32.  Good  Luck’ 
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